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Abstract 

We illustrate a completely analytic approach to Mel'nikov theory, 
which is based on a suitable extension of a classical method, and which 
is parallel and - at least in part - complementary to the standard pro- 
cedure. This approach can be also applied to some "degenerate" sit- 
uations, as to the case of nonhyperbolic unstable points, or of critical 
points located at the infinity (thus giving rise to unbounded orbits, 
e.g. the Keplerian parabolic orbits), and it is naturally "compati- 
ble" with the presence of general symmetry properties of the problem. 
These peculiarities may clearly make this approach of great interest 
in celestial mechanics, as shown by some classical examples. 



*Email: cicogna@df.unipi.it 
T Email: msantopr@math.uvic.ca 



1 Introduction 



It is certainly impossible to give a fairly complete list of the papers devoted 
to the applications of the Mel'nikov method ^], ||] for evaluating the onset 
of chaos arising in perturbed homoclinic (or heteroclinic) orbits. We will 
quote only some of the papers which are more directly connected with the 
present approach. 

The main purpose of this paper is to illustrate a completely analytic pro- 
cedure, based on a suitable extension of a classical method ||, which is 
parallel and - at least in part - complementary to the standard approach 
(see e.g. 0, [3], ^[). This procedure can be also applied to some "degener- 
ate" situations, as to the case of unstable nonhyperbolic points, or of critical 
points located at the infinity (thus giving rise to unbounded orbits, e.g. the 
Keplerian parabolic orbits), and it is naturally "compatible" with the pres- 
ence of general symmetry properties of the problem. For these reasons, apart 
from a clear "unifying" aspect, this method could be of great interest in ce- 
lestial mechanics, and it could be a contribution to the study of some of the 
questions about the onset of chaos in unbounded phase space systems and in 
the presence of unbounded orbits. It has been remarked indeed that in this 
situation chaos manifests itself in a particularly dramatic way J7| . 

Let us remark immediately that, in the above mentioned degenerate cases, 
i.e. in the lack of the hypothesis of hyperbolicity, standard results of pertur- 
bation theory cannot be directly applied; for instance, to preserve the critical- 
ity, we will have to impose a sufficiently rapid vanishing of the perturbation 
at the critical point. We can then extend the introduction of Mel'nikov func- 
tions, and show not only the existence of smooth solutions of the perturbed 
problem, approaching the critical points and playing in this context the role 
of stable and unstable manifolds, but also the possible presence of infinitely 
many intersections of these asymptotic sets on the Poincare sections, thus 
leading to a complicate dynamics typical of the homoclinic chaotic behaviour 



Being mainly interested in the methodological aspects, we will not devote 
special emphasis on new applications, but rather we will show, in the two 
last sections of this paper, how the method can be concretely applied in 
some typical situations, arising especially in celestial mechanics and general 
relativity. 



2 Statement of the Problem 

Although the procedure is completely general (indeed, we will state some 
results in a quite general setting in sections 3 and 4), we have actually in 
mind applications to celestial mechanics or to gravitational problems in gen- 
eral relativity; therefore we will restrict our attention mainly to problems 
described by Hamiltonians of the classical form 



where V(r) is some "long range" potential. Concretely, we have in mind 
potentials of the form 

„ (r) = _I^_<|_<| (2) 

where /3j are given parameters, which model several interesting situations, 
including standard Kepler problem, Kepler problems plus quadrupolar effects 
J7J, Manev problem ]ITJ, or the motion of a charged particle in the field of 



a Schwarzschild black hole in general relativity [|TT| [TI|, and many other 
situations. 

Let us write the equations of the motion, assuming for a moment that it 
lies in the plane 6 = ir/2, 

dV L 2 L 
r = p r , p r = - — + — - , ip = — , p v = (3) 
ar r J r z 

where we have put L = p^ = const ^ 0. We are interested in this paper in 
the appearance of chaotic behaviour related to the presence of homoclinic 
(or heteroclinic) orbits subjected to perturbations; therefore, the relevant 
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situations which may occur, depending on the values of the parameters $ in 
the potential (0), and which we are going to deal with in this paper, are the 
following: 

i) the presence of an unstable equilibrium point r u ^ for the first two 
equations in the system (|3D, which involve the variables r,p r . This point 
corresponds in the plane 9 = tt/2 to a unstable circular orbit 7 of radius r u . 
In this case, we have also the presence of a 1— parameter family of homoclinic 
bounded orbits biasymptotic to 7 (but see also case in) below). 

ii) the degenerate situation where the unstable equilibrium point is located 
at the infinity, i.e. r = 00, r = 0; the homoclinic orbits are in this case a 
family of parabolas. 

Hi) a "critical" case, with an unstable equilibrium point r u 7^ and another 
unstable equilibrium point located at r = 00, and a family of heteroclinic 
orbits connecting these points. 

It is clear that - due to the spherical symmetry of the problem - all 
conclusions and properties stated for the plane 9 = ir/2 are equally true for 
any plane for the origin in R 3 . 

Let us now choose and single out the following homoclinic (heteroclinic in 
case Hi)) orbit in the plane 9 = ti/2, denoted by x{t)i written in the spherical 
variables 

u:= (r,p r ,(p,p v ,6,p ) 

as follows 

%(t):= (R(t),R(t),$(t),L,n/2,0) (4) 

where R(t) and $(£) solve @ with the conditions i?(±oo) = r u in case i) or 
respectively R(±oo) = 00 in case ii), and R(— 00) = r u , R(+oo) = +00 in 
case in), and with $(0) = it. It is not necessary, for our purposes, to know 
explicitly the expression of the functions R(t) and $(£); it will be useful only 
to know that choosing R(0) = r (the turning point) in cases i) and ii), R(t) 
is an even function and $(£) an odd function of the time t. Let us remark 
that any other homoclinic orbit can be transformed by means of a rotation 
into the x(t) given by (f|). 

We now introduce a smooth (analytic) perturbation depending in general 
on all the variables u and time-periodic; the equations of the motion we are 
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considering are then, in general, 



p r + eg r (u,t) 

^2 n 2 



ri-i Pip d>V . , 

Pr = —+ 3 .2/3 _ T + e 9 P r[U, t) 

r 6 sm v dr 



$ = 2 PiP 2n +^(u,t) (5) 

r z sin 
= ^ + eg e {u,t) 



p"^ cos 9 
r 2 sin 3 



where e < 1. Let us also write (|5|) in a more convenient compact form 

u = f(u) + eg(u,t) (6) 

where 

/ = JV U H (7) 

and J is the standard symplectic matrix. Given any homoclinic orbit x(t) 
of the unperturbed problem, one has that \t '■= x{p ~~ satisfies, for all 
to £ R, 

^ = /<*.) (8) 

In order to find conditions ensuring the occurrence of intersections of 
stable and unstable manifolds of the critical point for the perturbed prob- 
lem, and hence the appearance of chaotic behaviour, we follow a (suitably 
extended) procedure which has been first used in this context (to the best 
of our knowledge) by Chow, Hale and Mallet-Paret Q in a problem with 
1 degree of freedom. We have first to look for smooth solutions near the 
homoclinic orbit; we then put 

u(t) = X (t - t ) + v(t - t ) (9) 

Substituting in @, we find for v(t) the equation (with the time shift t — > 
t + t ) 

v = A(t)v + eg(x(t),t + t ) + higher order terms in v and e (10) 
:= A(t)v + G(t,t ,e) 
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having used the shorthand notation G(t,t ,e), and where 



A(t):=(V„/)(x(t)) (11) 

This equation, or more often its first-order approximation 

v = A(t)v + eg(x(t),t + to) (12) 

is called the variational equation associated to (§) and to the orbit 

All the solutions v(t) of (TO) can be written, in implicit form (see |5|, [T3] ) 

v(t) = v h (t) + /' ^(s) G(s, t , e) cfe (13) 

where t\ is arbitrary, Vh(t) is any solution of the homogeneous linear problem 

v h = A(t)v h (14) 

and ^ is a fundamental matrix of solutions of ([14]) . We now have to look for 
the solutions of the homogeneous equation flT4]) and to control their behaviour 



for t — > ±oo. As well known [|], ||, [13| , this equation admits some solutions 
which remain bounded for all t G R and other solutions which diverge for 
t — > ±oo. The asymptotic behaviour is exponential in the case of hyperbolic 
unstable points, is like some power {t^ in the case of unstable point at r = 
oo, r = 0, and more in general of nonhyperbolic critical points. 

An obvious bounded solution is given by dx/dt. To find other solutions, 
the following results, which we will state for convenience in a quite general 
form, will be useful. 



3 General Results. The Role of Symmetry 

Let us start with the following definition (cf. ]T3, |T5[ ) 
Definition 1 A vector field 



S:=a(u)-V = ][>( M )A (15) 

i=l au i 
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is said to be a Lie-point symmetry (or - more exactly - the Lie generator 
of a Lie-point symmetry) for a general dynamical system (not necessarily 
Hamiltonian ) 

u = f(u) (16) 

if 

[<7.V,/-V] = (17) 
where [•, •] denotes the usual Lie commutator. 

We then can state: 

Lemma 1 Let x{t) be any given homoclinic orbit, solution of the problem 
(TBQ, and assume that this problem admits a Lie-point symmetry S = ct(m)-V. 
Then 

C(«) :=*(*(*)) 

is a solution of the homogeneous part of the variational equation asso- 
ciated to this problem. 

Proof. From direct calculation, using ( |TTD and (O), 

C = Jf(x(t)) = V«o- • x = (f ■ V„<7)(x(*)) = (V u / • a)( x (t)) = A(t)C. 



Remark. The obvious solution d\/dt of (|Tj) may be included into the 
solutions given by the above Lemma, indeed / • V = djdt is a (trivial) 
Lie-point symmetry of any autonomous dynamical system (Pp, expressing 
simply its time-flow invariance: if Uo(t) is a solution, then uo(t + to) is also a 
solution. 



Lemma 2 Given an Hamiltonian H , let K be a constant of motion for the 
Hamiltonian, i.e. {K,H} = 0. Then 

S:=JVK-V (18) 

is a Lie-point symmetry for the dynamical system u = JV U H . 

The proof is a straightforward verification. 
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4 General Results. 

Mel'nikov-type Conditions 

Still considering, from a general point of view, eq. as arising from an 
Hamiltonian with n degrees of freedom, assume that there are n bounded 
solutions C^'(i) an d n divergent solutions ry- a \t) of the homogeneous part 
(|T4| ) of the variational equation. Let us construct the fundamental matrix \I/ 
in ( |l3l) putting the solutions Q a * in the first n columns and the 77^ in the 
remaining columns of recalling that det \& = 1 and observing that \I/ 1 is 
a matrix having the 77^ in the first rows and — in the remaining rows, 
one easily sees (cf. ||) that ( |13D can be written (still in implicit form, see 

dD) 

v(t)=v h (t) + Ea ( [a] (t)£( V M(s)-JG(s,t ,e))ds (19) 
+ V [a] (t)j^(C [a] (s)-JG(s,t ,e))ds 

where Vh(t) is a linear combination of the C (*)> r^it). 

We now have to look for the existence of bounded solutions of (|10D ; more 
precisely we have to look for solutions v^~'(t) (and resp. v^ + '(t)) with the 
property of being bounded for t — > —00 (resp. t — > +00); these will provide 
precisely those solutions 

u^(t) = X (t-t )+v^(t-t ) (20) 

of @ which belong, by definition, to the unstable (resp. stable) manifold of 
the critical point. 

If we linearize (|T9| ) around the solution v(t) =0 (which amounts to delet- 
ing the higher-order terms in (|TUD, or to considering the variational equation 
in the form flT2]) ), and take into account the different asymptotic behaviour 
of the solutions C^(i) and rj^(t), it can be seen, as a consequence of the 
implicit-function theorem that the existence of bounded solutions 

both for t — > — 00 and t — > +00 is ensured if the following Mel'nikov-type 
conditions are verified 

M [a] (t )= [ + °°k [a] (t)- Jg(x(t),t + t ))dt = (21) 

J —00 v ' 
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Let us remark that if the perturbation is Hamiltonian, g = J^uW, as often 
happens in celestial mechanics, choosing C^CO = dx/dt, then the first of 
these conditions becomes 

M [1] (t )= f + °°{H,W}(x(t),t + t )dt = (22) 

J — oo 

Similarly, if, e.g., C' 2 ^) comes from a constant of motion K, according to 
Sect. 3, then 

M®(t )= / {K,W}( X (t),t + t )dt = (23) 

J — oo 

These conditions are identical to the conditions given e.g. in || ||, |I6| , where 
they are obtained by means of different procedures and hypotheses. It can be 
significant to remark that the present method then provides an extension of 
these formulas also to "degenerate" cases (nonhyperbolic points and possibly 
unbounded orbits), and to symmetries of more general nature, as in (p!5|). 

We have only to notice that, whereas the convergence of Mel'nikov inte- 
grals flSID is granted in the case of standard hyperbolic and isolated unstable 
points, the convergence is only "conditional", i.e. along a suitable sequence 
of intervals (see H for any details), in cases i) and in) of our classification 
in Sect. 2, and finally it must be controlled "by hand" in the non-hyperbolic 
case or in the case of critical point at the infinity. In the last cases one has to 
impose a sufficiently rapid vanishing of the contribution of the perturbation 
g(x(t),t) when t — > ±oo, i.e., as expected, when approaching the critical 
point. The precise rate of this vanishing will depend on the specific prob- 
lem in consideration (see Examples in the next sections, and also |17| for 



examples in 1 degree of freedom). For what concerns the regularity of the 
solutions and of the asymptotic manifolds, see e.g. [§, [T8|, [TOJ. 

Changing now the point of view, and considering the Poincare sections of 
the and solutions, the above arguments show that, once conditions 



(|2l| ) are satisfied, there occurs a crossing of the negatively and positively 
asymptotic sets on the Poincare section 0, |3], |(J . One usually imposes that 
the intersection is transversal; actually this condition is not strictly necessary, 
indeed it can be shown that it is sufficient that the crossing is "topological" 
30], i.e., roughly, that there is really a "crossing", from one side to the other, 
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but we do not insist on this point, which goes beyond the scope of the present 
paper. 

Thanks to the periodicity of the perturbation, one immediately deduces 
0, ^, f|, M that there is an infinite sequence of intersections of the positively 
and negatively asymptotic sets of the critical point in the Poincare section, 
leading to a situation typical of the homoclinic chaos. The presence of such 
infinitely many intersections is clearly reminiscent of the chaotic behaviour 
expressed by the Birkhoff-Smale theorem in terms of the equivalence to the 
symbolic dynamics of Smale horseshoes. Actually, this theorem cannot be 
directly used in the present context because its standard proof is intrinsically 
based on hyperbolicity properties [| . However, several arguments can be 
invoked even in the "degenerate" cases, which allow us to conclude that, if 
the conditions ( |2T| ) are satisfied, the perturbed problem exhibits a chaotic 
behaviour. We can refer e.g. to the classical arguments used in f2l| , and 
reconsidered by many others (see e.g. [0, [22|, ^3|, possibly resorting to 
singular coordinate transformations, such as the McGehee transformation in 



the case of critical point at the infinity, or the "blowing-up" method |25|, 
p6[ . More specifically, an equivalence to a "nonhyperbolic horseshoe" has 
been proved in ||, in which the contracting and expanding actions are not 
exponential but "polynomial" in time. The presence of Smale horseshoes 
and of a positive topological entropy has been also proved by means of a 
quite general geometrical or "topological" procedure |2(J which holds, in the 
presence of area-preserving perturbations, even in the nonhyperbolic case. 



5 Applications to Celestial Mechanics. 
Two Degrees of Freedom 

Coming back to the initial problems as stated in Sect. 2, the first step, ac- 
cording to the above discussion, is to write down the appropriate Mel'nikov 
conditions. It is clear that, apart from a rotation (this will change the expres- 
sion of the perturbation g in equations (f)]||) hito some new g which will also 
depend in general on the parameters of the rotation: this will be explained 
in detail later), we can always assume that the homoclinic orbit lies in the 
plane z = and in particular is given precisely by the orbit x(t) defined in 
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(^), Sect. 2. This greatly simplifies the variational equation, and in particular 
its homogeneous part (0); indeed it is immediate to verify that, with this 
choice, the two equations for the variations v 5, vq of 9 and pg are separated 
from the first four equations and admit regular bounded solutions; therefore 
the problem turns out to be 4-dimensional. An obvious symmetry for it is 
given by the rotations around the z— axis, generated by 

s = h (24) 

Then, as discussed in Sect. 3, two bounded solutions of the homogeneous 
part of the variational equation are 

C [1] = ^ = (R(t),R(t),m,o,o,o) (25) 

C [2] = (0,0,1,0,0,0) 

where R(t), $(t) have been defined in Sect. 2 (see (§,fD). We then get from 
([21] ) the two Mel'nikov conditions 

AfW = [ + ™\R(t)g P Mt),t + t )-R(t)g r (---) (26) 
J —00 L 

+ $(*)&„(...)]<** = 

/+00 
gpvixfot + h) dt = (27) 
-00 

Notice that usually one has g r = 0. In the case where the problem, included 
the perturbation, is completely planar, and that the perturbation is generated 
by an Hamiltonian W of the form 

W = W(r,<p,t) (28) 

then, under rotation, W is changed simply into W(r, ip + ip , t) and the above 
conditions become 



dr dip 



dt = 

(29) 
and 

«,) = / + °° ^(«M.«W+^' + «o) a = (30) 
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and it is easily seen that these can be also written according to the general 
form as given in ( f2~2] , f2~3| ) . 

As a first simple application of this situation, consider the classical Ke- 
pler problem with unstable equilibrium point at r = oo and r = and a 
perturbation not depending on <p, e.g. 

sin 2-Kvt 



W 

as in the classical Gylden problem 
ically satisfied, and the integral in 
function R(t j) 



(5 > 1/2) 



(31) 



|23|] . Then, condition ( p0|) is automat- 
) becomes (thanks to the parity of the 



cos2-Kut, 



R(t) 



sin27rz4 dt 



3 R(t) 8+1 

which clearly admits simple zeroes (the integral is easily seen to be 7^ 0). In 
this example, it is also simple to evaluate explicitly the asymptotic behaviour 
of the solutions of the variational equation: one has indeed M 



( = R(t) ~ \t 



-1/3 



rj{t) 



1 4/3 



for \t\ 



-00 



(32) 



and all conditions for the procedure given in Sect. 4 are satisfied, with - in 
particular - the condition 5 > 1/2 which ensures in this case the correct rate 
of vanishing of the perturbation at r = 00. Notice that this same condition 
would also guarantee that under the McGehee transformation [0, [ 
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the perturbation is not singular. Then chaotic behaviour is expected for this 
problem. 

A similar result holds for the more general (time-periodic) perturbations 
of the form (28) occurring e.g. in the restricted 3-body problems HT^ 



24 



These cases can be conveniently dealt with in this way. Assuming that 
W(R(t), $(£),£) -> for t -> ±00, then fl29|) (or (g|)) becomes 

-00 qw 



dt 



■(R(t),$(t) + ipo,t + t )dt = 



(33) 



(where clearly the derivative d/dt must be performed only with respect to the 
explicit time-dependence of W) ; introducing then the "Mel'nikov potential" 
W = W(y?o,£o) (cf- 123), corresponding to the perturbation W(r,cp,t): 



W(to,<Po):-- 



W(R(t),$(t) + <p ,t + t )dt 



(34) 
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then one gets from this definition and from 

M^(t , Vo ) = ^ = , M[ 2 ](t ,^ ) = ^ = (35) 

In other words, the two Mel'nikov conditions are equivalent to the existence 
of stationary points for the Mel'nikov potential W(<po,t ). On the other 
hand, W is a smooth doubly-periodic function, and such a function certainly 
possesses points to, Tfo where the two partial derivatives in fl35"D vanish, and 
this implies that the two conditions (|33|j30|) are certainly satisfied (one has 
only to check that these stationary points of W are isolated). 



6 Applications to Celestial Mechanics. 
Three Degrees of Freedom 

Let us consider finally a perturbation depending on both angles (p, 9. As 
already stated, to obtain Mel'nikov condition in the form (^6|,^j]), one has 
to transform the generic homoclinic orbit into the orbit x(t) given by (f|). 
This is obtained by means of a rotation defined by the following Euler angles 
(with the conventions and notations as in PS| ) : 



— Q, —i, —uj (36) 

where, with the language of celestial mechanics, i is the inclination of the 
plane of the orbit, u the angle of the perihelion with the line of nodes in the 
orbital plane, and Q is the longitude of the ascending node. 

Considering for simplicity a perturbation generated by an Hamiltonian 
not depending on the variables p, which we denote here, with a little abuse 
of notation, either by W(r, (p, 9, t) or by W(il, t), x = (x, y, z), and denoting 
by B the matrix of this rotation, the new expression W of the perturbation 
is obtained by replacing x with I?x. It is then easy to verify that Mel'nikov 
conditions (12611271) become 



/■+ 00 r • dW ■ /dW\ i 
MW( io , w ,Ji) = [ fl ^- + *(_)J A = (37) 

rm r +<xi /dW\ 

JK«fo, W ,n) = J_ x (^) t * = (38) 
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where 



/dW\ / \ 

(-— ) o = R(t) ( - G x {R{t), u, fi) sin + C 2 (. . .) cos $(*)) (39) 



<9y? /o 
with 

( dW\ 

Ci = (— — ) (cos a; cos — cos 2 sin a; sin fi) + 

f Q W\ , „ . . ON ( dW^ . . . 

+ — — (— cos u; sin Si — cos « sin a; cos i 2) + — — sin « sin u; 

,dW\ 

Co = ( t: — ) (sin u; cos Si + cos i cos sin Si) + 

( d W\ , . . n . ^ fdW^ . . 

+ — — (— sinwsinU + cos 2 cos u; cos S 2 ) — — — sinicosw 
\dx 2 Jo ^8x3^0 

and where (dW/dxi)o means that in the derivative of the given W with 
respect to Xi one has to replace x with I?x and finally put z = (or 9 = tt/2). 

If % = 0, i.e. if the problem is completely planar, including the perturba- 
tion, or if the perturbation is "generic" , i.e. has no "preferred" direction in 
the space (as often happens, see below for an example), and therefore it is 
not restrictive to choose i — 0, then one gets 

^ f dWs ,dWs . ^ ( 8W X . ,dWs 

Cl = (^)o COS ^ + (So Sm ^ ■ C - = (dx-)o Sm Vo-{^- 2 ) cosvo 

where (p = — (u + Si), and the above expressions (|^,^) are recovered. 

An example is provided by a problem in general relativity [II, [12j]. Con- 
sider indeed the motion of a relativistic charged particle in a gravitational 
field produced by a Schwarzschild black hole, perturbed by a homogeneous 
constant electric field. It can be shown [12] that the perturbation is given by 

W = F(r) (h cos ip sin 9 + l 2 sin ip sin 9 + l 3 cos 9) (40) 

where T{r) is suitable function, and 1 = (li,l2,h) is the direction of the 
perturbing field. This direction is generic and therefore no rotation is required 
(it would simply change the values of k, which are not fixed). 
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In this case, the first Mel'nikov condition in fl3"T|) is identically satisfied, 
because the perturbation is independent of the time (cf. (|33|) ), and the other 
one becomes 

/+oo 
J zr ( J R(t))cos$(t) dt (41) 
-oo 

which clearly admits simple zeroes, if only one can show that the integral is 
different from zero (see flTT||). 

Notice that in this example one has that the effective unperturbed po- 
tential admits bounded homoclinic orbits biasymptotic to an unstable point 
r u 7^ (which corresponds to unstable circles, as in case i) of our initial 
classification), therefore the above integral is expected to converge only con- 
ditionally §, 0. 

Similar results hold if the perturbation is produced by a constant homoge- 



neous magnetic field | I2| . It can be observed that, whereas the component of 
the electric and of the magnetic field on the plane of motion leads to a chaotic 
dynamics, the component normal to the plane does not. Since the problem 
is spherically symmetric, this argument can be applied to every plane for the 
origin. Thus, given an electric or magnetic field, on each plane for the origin 
(except at most the one normal to the field), chaos appears for a suitable 
choice of initial conditions. 
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